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Magnetoresistance in single crystals of black phosphorus is studied at ambient and hydrostatic 
pressures. In the semiconducting states at pressures below 0.71 GPa, the magnetoresistance 
shows periodic oscillations, which can be ascribed to the magneto-phonon resonance that is 
characteristic of high mobility semiconductors. In the metallic state above 1.64 GPa, the both 
transverse and longitudinal magnetoresistance show titanic increase with exhibiting super¬ 
posed Shubnikov-de Haas oscillations. The observed small Fermi surfaces, high mobilities 
and light effective masses of carriers in semimetallic black phosphorus are comparable to 
those in the representative elemental semimetals of bismuth and graphite. 


Black phosphorus (BP) is one of the stable allotropes of phosphorus at room temperature 
and consists of puckered honeycomb layers of phosphorus atoms in the ac-plane of the or¬ 
thorhombic crystal. As shown in the inset of Fig. 1(a), the directions along the zigzag chain 
and its normal correspond to the a- and c-axis, respectively. The monolayer of BP (called 
phosphorene) is known as a high mobility semiconductor having a direct gap of 2 eV at the F 
point and attracts recent attention as a new candidate for electrical devices next to graphene. 
In bulk BP, the honeycomb planes stack along the h-axis and are weakly bonded with the van 
der Waals force with alternatively shifting along the c-axis, i.e., forming ABAB- • • stacks. 
The interlayer interaction involves band dispersion along the F-Z line. As a result, the bulk 
BP becomes a narrow gap semiconductor with the direct gap of 300 meV at the Z point. 
Application of hydrostatic pressure to BP further increases the interlayer interaction and then 
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closes the gap at about 1.5 GPa. This pressure-indueed semiconductor-semimetal (SC-SM) 
transition has been experimentally observed in the transport and the infrared absorption mea¬ 
surements under pressure.Although the semimetallic state in this high mobility material is 
expected to exhibit unusual quantum transport properties, the details have not been revealed. 
This motivates us to earry out the magnetoresistance measurements in BP under pressure. 

Single erystals of BP were synthesized under high pressure.^^ Transverse and longitudi¬ 
nal magnetoresistance were measured by the four-probe method on samples with the typical 
dimensions of 1.5 x 1.0 x 0.1 mm^. Gold wires of 30 /rm-thiek were attached to the eleaved 
surfaees of the samples with MRX-713J earbon epoxy. High pressure environments up to 
2.4 GPa were realized by a piston eylinder-type pressure cell. The pressures in the sample 
space were determined by monitoring the supereonducting transition temperatures of Pb set 
together with the samples. Glycerin or Daphne 7373 were used as pressure media. For mag¬ 
netoresistance measurements at ambient pressure, we utilized non-destructive 56-T magnets 
at The Institute for Solid State Physics at The University of Tokyo with the duration of about 
36 ms, whereas the measurements under pressure were earried out up to 14 T with using su¬ 
pereonducting magnets of Physieal Property Measurement System (Quantum Design, Ine.). 
Magnetic fields are applied along the a-axis throughout the entire measurements. 

Figure 1(a) shows the temperature dependenee of the resistivity along the c-axis (pc) at 
various pressures. Above 1.34 GPa, pc shows metallie behavior in the whole temperature 
range between 2 and 300 K. First, we focus on the transport properties of the semiconducting 
BP under ambient and relatively low pressures. Figure 1(b) shows the field dependenee of pc 
at 77 K with a magnetic field (B) up to 56 T applied along the a-axis. The pc shows positive 
magnetoresistanee aceompanied by some hump struetures below 30 T. These struetures are 
periodic in terms of B~^ as shown by the solid and open eireles in Fig. 1(e). This periodic 
structure has been aseribed to the magneto-phonon resonanee (MPR).^^ The MPR is the res¬ 
onant seattering of earners by optical phonons and observed as enhaneement of resistanee 
when the following resonance eondition is fulfilled: 

nojo = moj^ (N = 1, 2, 3,...) (1) 

where tiioo is the energy of the optieal phonon and a>c = eBIm* is the cyelotron frequency. 
The e and m* represent elemental eharge and effeetive mass of earners, respeetively. The 
MPR beeomes hardly visible at low temperatures beeause the thermal excitation of the optical 
phonon is suppressed. Equation (1) ean be transformed to the relation between inversed field 
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Fig. 1. (Color online) (a) Temperature dependence of pc at various pressures from 0.0 to 2.37 GPa. The inset 
illustrates the relation between the crystal structure and applied currents, (b) Magnetic field dependence of pc at 
ambient pressure and temperature of 77 K. Magnetic field up to 56 T was applied along the a-axis. (c) Second 
derivative of pc as a function of at ambient and hydrostatic pressures. Peak and dip structures at each 
pressure are denoted by solid and open symbols, respectively, (d) Relation between inversed peak/dip fields B ' 
and the index N in Eq. (1) (see text). 


and the integer index as 


R = —N. (2) 

B rrfojQ 

With using coq = 2.60 x 10'^ s“' (138 cm“^), which is the energy of the longitudinal opti¬ 
cal (LO) phonon determined by the Raman scattering,m* estimated from Eq. (2) is 0.153 
mo at ambient pressure, which is consistent with the previously reported values of 0.146 mo 
at 20 and 0.166 mo at 30 determined by experiments of the cyclotron resonance. 
Observation of the MPR indicates the existence of well-defined Landau levels in BP, i.e., 
ojcT = juR > 1 above 4 T. Here, r and /r denote the relaxation time and mobility of carriers, 
respectively. Therefore, /i > 0.25 m^ s“^ = 2500 cm^ s“' at 77 K, which is consistent 
with the earlier report. Figure 1(c) shows the pressure dependence of the second deriva¬ 
tive -d^c/dR^ at 77 K. The period of oscillation becomes long with applying pressure, which 
is interpreted as the increase of the coefficient e/(m*mo) in Eq. (2). Assuming the relevant coo 
is insensitive to the pressure as the other phonon modes are,^®^ we can evaluate the pressure 
dependence of m* from the slopes of Fig. 1(d). In Fig. 1(d), half-integer indices correspond 
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Fig. 2. (Color online) (a) Transverse magnetoresistance at 1.64 GPa. (b) Second derivative of pc shown in (a). 
Solid and open arrowheads denote the oscillatory structures originate from different cross-section of the Fermi 
surface, (c) Fast Fourier transform (FFT) of SdH oscillation shown in (b). Detected frequency components (a, 
P, and 2a) are denoted by solid arrows. The inset shows the temperature dependence of FFT magnitude at a and 
p. Dashed lines are fitted curves with using Eq. (3). (d) The temperature dependence of pc in several magnetic 
helds at 1.64 (solid lines) and 2.00 GPa (broken lines). 


to the local minima in Fig. 1(c). Figure 1(d) clearly indicates the reduction of m* by applied 
pressure. According to the A: • /? theory, the energy gap is roughly proportional to m* in the 
case of direct band gap semiconductors. Hence, this implies the suppression of the energy 
gap by pressure. 

Now, we present the detailed transport properties of semimetallic BP. Figure 2(a) illus¬ 
trates Pc at 1.64 GPa as a function of magnetic field along the a-axis. At 2 K, the magne¬ 
toresistance normalized by the value at zero-field reaches pc(B)/pc(0) ~ 1000 at 14 T. Such 
a drastic change in magnetoresistance has been known in bismuth^^^ and graphite,typical 
elemental semimetals, and attracts renewed interests in the other semimetals. 

In addition to this approximately linear magnetoresistance, we can identify the super¬ 
posed modulation. The modulated components are clearly visible in their second derivatives 
shown in Fig. 2(b). Contrary to the MPR in the semiconducting state, the oscillating compo¬ 
nent periodic to grows up as the temperature decreases, and hence, can be ascribed to the 
Shubnikov-de Haas (SdH) oscillations. We can see a major frequency component (solid ar¬ 
rowheads) and additional small structures (an open arrowhead) in Fig. 2(b), suggesting there 
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are at least two frequeney eomponents. Frequeney speetrum and its temperature dependenee 
are shown in Fig. 2(o). We ean see two distinet peaks, a (4.9 T) and >0(11 T), and the seeond 
harmonie 2a (9.8 T) as marked by arrows. Assuming the semimetallie state, there should be 
at least two Fermi poekets, i.e., eaeh one eleetron- and hole-poeket. The eross-seetional areas 
of eaeh Fermi surfaee are 0.47 xlO~^A“^ and 1.0 xlO~^A^^ for a and yS, whieh is approxi¬ 
mately 0.017 and 0.037 % of that of the first Brillouin zone at the = 0 plane. For spherieal 
Fermi surfaees, these small poekets roughly eorrespond to the earner densities of 10^^ em“^. 
Although the aetual Fermi surfaee should be highly anisotropie,^^’^^^ its earner density ean 
be eomparable to those of bismuth and graphite. 

In the semielassieal treatment, the SdH oseillations are represented by Lifshitz-Kosevieh 
(LK) formula.^®’^^^ In this formula, the thermal damping faetor of amplitude Rj is given as 
follows: 

^ _ In^kBTIfiojc 

^ sinh(2;r^kB7’/^<^c)’ 

where k^ is the Boltzmann eonstant. The observed temperature dependenee of the oseillation 
amplitude ean be reprodueed by Eq. (3) with m* ~ 0.02 mo for both a and yS as shown in the 
inset of Fig. 2(o). 

Next, we demonstrate the magnetoresistanee at various pressures in the semimetallie state. 
Figure 3(a) shows Pc at 2.0 K as a funetion of applied field along the a-axis. As illustrated 
in the FFT speetra [Fig. 3(b)], the peaks a and yS deteeted at 1.64 GPa move toward higher 
frequeney, indieating the enlargement of the Fermi surfaees by pressure. We investigated the 
pressure dependenee of these peaks in another sample as shown by the open symbols in Fig. 
3(e) and obtained eonsistent results. The pressure dependenee of m* determined through anal¬ 
yses of the MPR and the SdH is shown in Fig. 3(d). m* of semimetallie region is signifieantly 
smaller than that of semieondueting region. 

At the vieinity of the SC-SM transition pressure, the quantum limit state in whieh all 
the earriers are aeeommodated in the lowest Landau subband ean be realized at around 10 T. 
Reeent experimental studies on graphite suggested the emergenee of an exeitonie phase at the 
vieinity of the quantum limit state in high magnetie fields.To elarify whether an exotie 
phase emerges in BP, we studied the temperature dependenee of the resistivity as shown in 
Fig. 2(d). The temperature dependenee at zero-field is metallie as mentioned above, while 
in magnetie fields, they turn to the semieondueting ones in temperature up to 200 K. Sueh 
apparent semieondueting-like behavior in semimetals ean be explained by simple two earner 
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Fig. 3. (Color online) (a) Magnetic field dependence of pc at various pressures, (b) FFT spectrum of the SdH 
oscillations under several pressures. The a and j3 peaks are traced by solid and open arrowheads, respectively, 
(c) Pressure dependence of the a and /? peaks. Open symbols indicate the result on another sample piece, (d) 
Pressure dependence of the effective mass estimated from MPR in semiconducting region (<1.5 GPa) and from 
the temperature dependence of SdH amplitude in semimetallic region ( > 1.5 GPa). Vertical errorbars represent 
the statistical error in curve fitting. 


model.In the semielassieal approaeh, the transverse resistivity (p±), is expressed as^®^ 

_ PePh(Pe + Ph) + (Pe-^h + Ph^e)^^ 

“ (Pe+Ph)" + (i?e+i?h)252 ’ 

where i?e,h = 1 / («e,h<?e,h) is the Hall eoeffieients of eleetron and hole earners (ne.h and (?e,h 
represent the density and the eharge of eleetron and hole earriers, respeetively) and pe,h is the 
resistivity of eaeh earrier at zero-field. In a eompensated semimetal, we ean set -qg = and 
Ug = Hh, i.e., -Rg = Ri^= R. Assuming pg(T) = ph(T) = po(T) at zero-field, Eq. (4) is reduced 
to 


i2 d2 


PoiT) ^ R^B 


Px - —^5— + —777- (5) 

2 2po(T) 

In low-carrier systems, the second term in Eq. (5) becomes dominant, and hence px can show 
the superficial insulating-like behavior even in a metallic sample. 

To clarify the intrinsic features of Pg{T) or p\i{T), it is rather suitable to study longitudinal 
magnetoresistance. Eigure 4(a) shows a-axis resistivity (pa) as a function of B applied parallel 
to the a-axis at pressure of 1.65 GPa. In this longitudinal geometry, we clearly observed 
SdH oscillations with the frequencies almost identical to those in the transverse one at 1.64 
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Fig. 4. (Color online) (a) SdH oscillation in the longitudinal resistance pa at 1.65 GPa. (b) Temperature de¬ 
pendence of Pa at 1.65 GPa under several magnetic fields. Horizontal axis is common to (c). (c) Temperature 
dependence of pa normalized by the resistivity at zero field pa(0) at 1.65 GPa. Fitting curves based on Eq. (6) 
in the text are shown by broken lines, (d) Fan diagram at a peak for transverse (B ± I) and longitudinal (B || I) 
configurations. 


GPa. Thereby, quantum limit state is aehieved at fields above 11 T. Even in this longitudinal 
geometry, we observed non-saturating signifieant positive magnetoresistanee. At 2 K, the 
ratio Pa(5)/pa(0) reaehes approximately 50 at 14 T. 

As shown in Fig. 4(b), this longitudinal magnetoresistanee shows peak strueture at around 
30 K as a funetion of temperature. Sueh non-monotonie behavior, however, may not be as- 
eribed to the emergenee of a novel phase below this temperature. Figure 4(e) shows the tem¬ 
perature dependenee of the Pa(5) normalized by the value at zero field Pa(0). The traees show 
monotonie behaviors at all the fields. Aeeording to a transport theory in the quantum limit,^^^ 
The longitudinal resistivity p\\{B, T)lp\\(0, T) for the elassieal statisties, in whieh the Fermi 
energy is smaller than represented as, 

P\\{B, T) _l I hcoc \ _ fie B 
PlliO,T) ~ 3\kBTl~ 3k^m*T 

in ease that the seattering meehanism is dominated by d-funetion impurity potential or aeous- 
tie phonon. As shown with broken lines in Fig. 4(e) , Pa(5, T)lp^(0, T) is reasonably repro- 
dueed by Eq. (6) at all the field by adjusting single parameter m* = 0.005 mo above 20 K 
where the elassieal statisties likely being relevant. Although the m* value used in this analysis 
is different from that estimated above, this diserepaney is not unusual with taking into aeeount 
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the crudeness of the model. Therefore, unusual temperature dependence of the transverse and 
longitudinal magnetoresistance in semimetallic BP does not show the emergence of a gapped 
phase at high magnetic field, but can be simply understood as characteristic behavior in a 
high mobility and low carrier-density semimetal comparable to bismuth and graphite. 

Here, let us make a comment on the phase of the SdH oscillations. In Fig. 4(d), we plot 
the relation between inversed dip fields of Pa and pc and the Landau index. Similar analyses 
are frequently utilized in topological materials to evaluate the Berry phases from the values 
of the horizontal intercept in this diagram. The intercept for pc is 0.00 + 0.01 similar to 
two-dimensional Dirac system, whereas that for pa is 0.15 + 0.02 even for the same field 
direction. In addition, we note that such analyses on resistivity make sense only in the case 
that Hall conductivity is much larger than the diagonal one.^*^ Therefore, we cannot make 
reliable argument from the analyses of this diagram. 

Finally, we show the results of band calculation under pressure and compare them with 
the experimental results. Figure 5 (a) shows the band structure at pressure of 2.0 GPa by the 
first-principle calculations within the generalized gradient approximation using the OpenMX 
code.^°^ Two and one primitive obitals for s- and p-obritals are used for the pseudo-atomic 
orbitals basis functions, respectively. For the crystal structure parameters at 2.0 GPa, we lin¬ 
early extrapolated the values at ambient pressure and 0.8 GPa of Ref. 4. (We have checked 
that we can obtain almost the same results as Fig. 5 (a) even by using the crystal structure 
parameters given in Ref. 32.) 

At ambient pressure, the obtained band structure (not shown) is consistent with previous 
works^’^^^ and the direct band gap at the Z point is 130 meV. At pressure of 2.0 GPa, the band 
overlap between the conduction and valance bands at the Z point is 140 meV. There are one 
hole pocket at around the Z point and four electron pockets locating between the F-A line. 
(Their energy minimum are slightly off the F-A line.) The hole pocket is highly anisotropic: 
the ratio of Fermi wavevectors is : k^~^ = 21:13:3, while the electron pockets 

are rather isotropic. Although the volume of the hole pocket should be four times larger than 
that of the electron pockets, the cross-sectional area of hole pocket is not necessarily so due 
to the highly anisotropic hole pocket. (The a-axis corresponds to the Z-A direction.) 

At the final stage in preparing this manuscript, we became aware of the paper has been 
submitted to arXiv that reports similar quantum transport in semimetallic BP.^^^ 

In conclusion, we investigated the magnetoresistance in semiconducting and semimetallic 
black phosphorus under hydrostatic pressures. In the semiconducting state, we observed the 
magneto-phonon resonance that is characteristic of high mobility materials. In the semimet- 
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Fig. 5. (Color online) (a) Band structure at pressure of 2.0 GPa. (b) Brillouin zone and some of the symmetric 
points, (c) Fermi surfaces: one hole pocket at around the Z point and four electron pockets between the F-A line. 
The b-3xis corresponds to the F-Z direction and the c-axis does to F-A direction. 


allc State above 1.6 GPa, we observed clear Shubnikov-de Haas oscillations and positive giant 
magnetoresistance effects. The analysis of the oscillation revealed its extremely small carrier 
density and light effective mass, which suggests that black phosphorus can be a novel material 
to investigate the exotic quantum physics expected in the quantum limit state. 
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